Resonant activation in bistable semiconductor lasers 
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We theoretically investigate the possibility of observing resonant activation in the hopping dynam- 
ics of two-mode semiconductor lasers. We present a series of simulations of a rate-equations model 
under random and periodic modulation of the bias current. In both cases, for an optimal choice of 
the modulation time-scale, the hopping times between the stable lasing modes attain a minimum. 
The simulation data are understood by means of an effective one-dimensional Langevin equation 
with multiplicative fluctuations. Our conclusions apply to both Edge Emitting and Vertical Cavity 
Lasers, thus opening the way to several experimental tests in such optical systems. 
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I. INTRODUCTION 

It is currently established that stochastic fluctuations 
may have a constructive role in enhancing the response 
of nonlinear systems to an external coherent stimulus. 
Relevant examples are the enhancement of the decay 
time from a metastable state (noise-enhanced stability) 
0, 0| , the synchronization with a weak periodic input sig- 
nal (stochastic resonance) Q or the regularizaton of the 
response at an optimal noise intensity (coherence reso- 
nance) [Hi- 

Another instance is the phenomenon of resonant acti- 
vation that was discovered by Doering and Gadoua ||. 
They showed that the escape of an overdamped Brown- 
ian particle over a fluctuating barrier can be enhanced 
by suitably choosing the correlation time of barrier fluc- 
tuations themselves. In other words, the escape time 
from the potential well attains a minimum for an optimal 
choice of such correlation time. Since its discovery, the 
phenomenon received a considerable attention from the- 
orists (see e.g. Refs. H, 0, S, S E3|)- Detailed studies by 
means of analog simulations have also been reported for 
both Gaussian and dichotomous fluctuations [ll|. More 
recently the phenomenon has been shown to occur also 
for the case in which the barrier oscillates periodically 



To our knowledge, experimental evidences of resonant 
activation were only given for a bistable electronic cir- 
cuit [HI and, very recently, for a colloidal particle sub- 
ject to a periodically-modulated optical potential [lo] ]. 
It is therefore important to look for other setups where 
the effect could be studied in detail. As a matter of 
fact, multimode laser systems are good candidates to 
investigate noise-activated dynamics like the switching 
among modes induced by quantum fluctuations (sponta- 
neous emission) [16 j . In particular, semiconductor lasers 
proved to be particularly versatile for detailed experi- 



mental investigations of modulation and noise-induced 
phenomena like stochastic resonance [13, EH and noise- 
induced phase synchronization [l9| . In those previous 
studies, the resonance regimes are attained by a suit- 
able random modulation of the bias current which can 
be tuned in a well-controlled way. It is thus natural to 
argue about the possibility of observing resonant activa- 
tion with the same type of experimental setup. 

In this paper, we theoretically demonstrate the phe- 
nomenon of resonant activation in a generic rate- 
equations model for a two-mode semiconductor laser un- 
der modulation of the bias current. The basic ingredients 
that act in the theoretical descriptions are a fluctuating 
potential barrier and some activating noise. In the laser 
system, the latter is basically provided by spontaneous 
emission while current fluctuations, that appear addi- 
tively into the rate equations, effectively act multiplica- 
tively if a suitable separation of time scales holds [2(|. In 
a previous paper [2l| . we have explicitely demonstrated 
such multiplicative-noise effects on the mode-hopping 
dynamics. This was shown by a reduction to a bistable 
one-dimensional potential system with both multiplica- 
tive and additive stochastic forces. Several predictions 
drawn from such a simplified model are in good agree- 
ment with the experimental observations carried out for 
a bulk, Edge-Emitting Laser (EEL) [2l| . In the present 
context, we will show that this reduced description is of 
great help in the interpretation of simulation data. 

The outline of the paper is the following. In Sec. II 
we recall the model for a two-mode semiconductor laser. 
In Sec. Ill we present the numerical simulation for two 
physically distinct cases displaying resonant activation. 
These results are discussed and interpreted by compar- 
ing with the reduced one-dimensional Langevin model 
mentioned above (Sec. IV). We draw our conclusions in 
Sec. V. 



II. RATE EQUATIONS 
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Our starting point is a stochastic rate-equation model 
for a semiconductor laser that may operate in two longi- 
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tudinal modes whose complex amplitudes are denoted by 
E± . Both of them interact with a single carrier density 
iV that provides the necessary amplification. The two 
modes have very similar linear gains, provided that their 
wavelengths are almost equal and they are close to the 
gain peak. Let J(t) denote the bias (injection) current, 
the model can be written as \2m 



(1 + ia)gj, 



(1 + ia)g. 



1 



E+ + y/2D sp Nb 



£L + y/2D ap N£- 



N = 7 



J(t) ~ N - g + \E + \ 2 - g_\E^ 



(la) 



(lb) 



(lc) 



where 7 is carrier density relaxation rate, a is the 
linewidth enhancement factor [22j • The modal gains read 



9± 



N ± e(N - N c ) 
l + s\E±\< 



(2) 



where e determines the difference in differential gain 
among the two modes while iV c defines the carrier 
density where the unsaturated modal gains are equal. 
The parameters s and c are respectively the self- and 
cross-saturation coefficients. The £± are two indepen- 
dent, complex white noise processes with zero mean 
[<&=(*)) = 0] and unit variance = M(t-f')] 

that model spontaneous emission. The noise terms in 
Eqs. (fTaj) and (jlb[) are gauged by the spontaneous emis- 
sion coefficient D sp . 

All quantities are expressed in suitable dimensionless 
units. In particular, time is normalized to the photons' 
lifetime, which for semiconductor laser is typically of the 
order of a few picoseconds or less (see e.g. [22|, [23, [24]]) 

A detailed analysis of the stationary solutions of 
Eqs. ((T|) is reported in Ref. (25[. For a constant bias cur- 
rent J(t) = Jo and D sp = 0, Eqs. |T]) admit four different 
steady state solutions: the trivial one E± = 0, two single- 
mode solutions — E+ 7^ 0, E- — and viceversa — and 
a solution where both modes are lasing, E± 7^ 0. For 
N c > 1, and c > s, there exist a finite interval of Jo val- 
ues for which the two single-mode solutions coexist and 
are stable while the E± 7^ is unstable (bistable region) . 
Here, for D sp > the laser performs stochastic mode- 
hopping, with the total emitted intensity remains almost 
constant while each mode switches on and off alternately 
at random times. We point out that the emission in each 
mode is nonvanishing even in the "off" state, as the av- 
erage power spontaneously emitted in each mode at any 
time is given by 4D sp N [recall that Eqs. (H|) are usually 
interpreted in Ito sense [HI]. Observation of this be- 
haviour has been rep orted in several experimental works 
on EELs [II, III, If. 

We remark that while Eqs. ([T]) aim at modeling EELs, 
the results presented henceforth would apply also to po- 
larization switching in Vertical Cavity Surface Emitting 



Lasers (VCSELs). Indeed, experimental data [2i| show 
strong similarities between this phenomenon and the lon- 
gitudinal mode dynamics. On the theoretical side, this 
analogy is supported by the fact that the polarization dy- 
namics in VCSELs is described by models that are math- 
ematically similar to the one discussed here [3(| HH, HH . 

In the following, we will focus on the effect of the ex- 
ternally imposed fluctuation/modulation of the injected 
current. This situation is modeled by letting 



J{t) = J +6J(t) 



(3) 



The DC value Jo sets the working point and will be al- 
ways chosen to be in the bistability region. We focus 
on the case in which SJ is a Ornstein-Uhlenbeck process 
with zero average (SJ(t)) = and correlation time r: 



SJ 



SJ_ 

T 



2D, 



that means 



(SJ{t)SJ{0)} = Djexp(-\t\/T) 



(4) 



(5) 



This choice is suitable to model a finite-bandwith noise 
generator. Notice that r and the variance of fluctuations 
Dj = (SJ 2 ) can be fixed independentely. 

Another case of experimental interest that we will con- 
sider is using the current modulation 



SJ 



A sin Qt 



(6) 



To assess the nature of the stochastic process at hand, 
it is important to introduce the relevant time scales. We 
define first of all the switching or relaxation time Tr as 
the typical time for the emission to change from one mode 
to the other. The main quantities we are interested in are 
the Kramers or residence times T± defined as the average 
times for which the emission occurs in each mode. In 
semiconductor lasers T± are generally much larger than 
Tr. Typically, Tr ~ 1 — 10ns while residence times may 
range between 0.1 and 100 /is [H, |33|. The third time- 
scale is of course given by the characteristic time of the 
external driving, namely, r and 2tt/Q respectively. 

In the following, we will study how the hopping dy- 
namics changes upon varying these latter parameters as 
well as the strength of the perturbation. 



III. NUMERICAL SIMULATIONS 

In this Section we present the outcomes of a series of 
numerical simulation of Eqs. ([1]). In Ref. [2l| it was 
observed that the sensitivity of each of the T± on the 
imposed current fluctuations may be notably different 
depending on the parameters' choice. This is a typi- 
cal signature of the multiplicative nature of the stochas- 
tic process. In particular, one can argue [2~ll ] that such 
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"simmetry-breaking" effects mostly depend on the ratio 
serf 5 where 



The parameter a represents the gain saturation induced 
by the total power in the laser, while S describes the 
reduction in gain saturation due to partitioning of the 
power between the two modes. 

The possibility of obtaining qualitatively different re- 
sponses depending on the actual parameters corresponds 
to the different ex per imental observations reported for 
both EELs [II Hill and VCSELs [13, HI- Those two 
classes of lasers were indeed found to display markedly 
different simmetry-breaking effects under current modu- 
lation. To account for those features, we consider two 
different sets of phenomenological parameters. For defi- 
niteness, in both cases we fix e = 0.1, s = 1.0, N c = 1.1, 
7 = 0.01 and change the values of c and D sp (see Table 
I). The first set (6 — 0.05) corresponds to the case in 
which added modulation changes the hopping time scale 
in an almost symmetric way. On the contrary, in the 
second case (5 = 0.15) the asymmetry effect of the noise 
is stronger [2lJ. We can thus consider the two as rep- 
resentative of the VCSELs and EELs case respectively. 
The value of Jo has been empirically adjusted to yield 
T + ~ T_ = T s and an almost symmetric distribution of 
intensities in absence of modulation. The actual values 
are about 10% above the laser threshold. The sponta- 
neous emission coefficient D sp has been chosen to yield a 
value of the residence times of the same order of magni- 
tude of the experimental ones. 

In the following, we decide to set a = which is ap- 
propriate for our EEL model where the phase dynamics 
is not relevant [21[ . This choice may however not be fully 
justified for the VCSEL case. In this respect, the sim- 
ulations presented below are representative of the VC- 
SEL dynamics only in a qualitative sense. Nonetheless, 
it should be pointed out that a ID Langevin model in- 
dependent of a describes also the VCSEL case (30L l3lj|. 
Since resonant activation is mainly due to the multiplica- 
tive noise effect described by such equations [see Eq. J9|) 
below] we consider this as an indirect proof that phe- 
nomenology we will report below should be observable 
also in the VCSEL case. 

The largest part of the simulations were performed 
with Euler method with time steps 0.01-0.05 for times 
in the range 10 7 — 10 s time units depending on the val- 
ues of r and f2. For comparison, some checks with Heun 
method [34| have also been carried on. Within the sta- 
tistical accuracy, the results are found to be insensitive 
the the choice of the algorithm. 

A. Stochastic modulation 

Let us start illustrating the results in the case of 
stochastic current modulation (Eq. (01). In Fig. Q] and 



TABLE I: The parameter values used in the two series of 
simulations of Eqs. (JTJ, the other values are given in the text. 



c 


D S p 


Jo 


5 


a 


1.1 


0.7 x 10" 5 


1.197 


0.05 


1.05 


1.3 


1.5 x 10" 5 


1.194 


0.15 


1.15 



[2] we report the measured dependence of the residence 
times T± on the correlation time r for the two parameter 
sets given in Table I and different values of the noise vari- 
ance Dj. In all cases, the curves display well-pronounced 
minima at an optimal value of r. This is the typical sig- 
nature of resonant activation. The minima are almost 
located between the relaxation time Tr and the hopping 
time T s (marked by the vertical dashed lines). The val- 
ues of Tr reported in the figures have been estimated 
from the reduced model discussed in the next Section, 
see Eq. (fT4"|) below. 

The effect manifest in a different way for the second 
parameter set. In the case of Fig. [1] both times attain a 
minimum, albeit with different values. On the contrary 
the data of Fig.[2]show that one of the two times is hardly 
affected from the external perturbation regardless of the 
value of t. In other terms, we can tune the current corre- 
lation in such a way that emission along only one of the 
two modes is strongly reduced (about a factor 10 in the 
simulation discussed here). 



B. Periodic modulation 

Let us now turn to the case of sinusoidal current mod- 
ulation (Eq. [B]). In Fig. and 2] we report the measured 
dependence of the residence times T± on the frequency O 
for the two parameter sets given in Table I and different 
values of the amplitude A. For comparison with the pre- 
vious case we choose A such that the RMS value of (|B|) 
is roughly equal to the variance of i.e. A ~ y/2Dj. 

As in the previous case, the curves display resonant 
activation at an optimal value of f2. For the second set 
of parameters, one of the two hopping times is more re- 
duced than the other (compare Fig. 0] with Fig. It 
should be also noticed that the data in Fig. [2] display 
some statistical fluctuations while the curves for the pe- 
riodic modulation are smoother. 



IV. INSIGHTS FROM A REDUCED MODEL 

In order to better understand the activation phe- 
nomenon it is useful to reduce the five-dimensional dy- 
namical system ([1]) to an effective one-dimensional sys- 
tem. This has been accomplished in Ref. [2lfl . For com- 
pleteness, we only recall here some basic steps of the 
derivation. In the first place, we introduce the change of 
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FIG. 1: (Color online) Simulations of the rate equations with 
Ornstein-Uhlenbeck current fluctuations, parameter set with 
c = 1.1 (see text and Table I): residence times T+ (squares) 
and T- (circles) for increasing values of the current variance 
D.j. The values of the relaxation time Tr and the hopping 
time T a (in absence of modulation) are marked by the vertical 
dashed lines. 
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FIG. 2: (Color online) Simulations of the rate equations with 
Ornstein-Uhlenbeck current fluctuations, parameter set with 
c = 1.3 (see text and Table I): residence times T+ (squares) 
and T_ (circles) for increasing values of the current variance 
Dj. 



coordinates 

E + = r cos 4> exp itp + , E— = r sin exp iip_ . (8) 

In these new variables, r 2 is the total power emitted by 
the laser, and <f> determines how this power is partitioned 
among the two modes. The values (j> — 0, tt/2 correspond 
to pure emission in mode + and — respectively. The 
phases ip± do not influence the evolution of the modal 
amplitudes and carrier density and can be ignored. 

In order to simplify the analysis, we assume that (i) 
The difference between modal gains is very small, i.e., 
Ar c >l,£<Sl,c>s; (ii) the laser operates close enough 
to threshold, so that r 2 «l and the saturation term is 
small: in this limit, r and N decouple to leading order 
from </>; (iii) r and N can be adiabatically eliminated and 
(iv) only their fluctuations around the equilibrium val- 
ues due to J are retained. This last assumption holds for 
weak spontaneous noise and amounts to say that r and 



N are stochastic processes given by nonlinear transfor- 
mations of J (see Eqs. (16) in Ref. [U)- This requires 
that J does not change too fast. For example, in the 
case of the Orstein-Uhlenbeck process, Eq. t should 
be larger than the relaxation time of the total intensity. 
The validity of the above reduction has been carefully 
checked against simulations of the complete model [2l| . 
For the scope of the present work, we performed a fur- 
ther check by comparing the spectrum of fluctuations of 
r 2 with the imposed one, Eq. (H|). Indeed, the behaviour 
is the same for r > Tr while for shorter r some differ- 
ences are detected. This means that the reduced descrip- 
tion discussed below becomes less and less accurate. On 
the other hand, in this regime spontaneous fluctuation 
should dominate and this limitation become less relevant 
for our purposes. 

Altogether, the hopping dynamics is effectively one- 
dimensional and is described by the slow variable <p. Its 
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Modulation period 2n/Q. 



FIG. 3: (Color online) Simulations of the rate equations with 
sinusoidal modulation of the current, parameter set with c = 
1.1 (see text and Table I): residence times T+ (squares) and 
T_ (circles) for increasing values of modulation amplitude A. 
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FIG. 4: (Color online) Simulations of the rate equations with 
sinusoidal modulation of the current, parameter set with c = 
1.3 (see text and Table I): residence times T+ (squares) and 
T_ (circles) for increasing values of modulation amplitude A. 



evolution is ruled by the effective Langevin equation 

2D, 



a cos 2(f) + b 



sin 2(j) 



tan 2(f) 



y/ZDtt* (9) 



where, together with ([7]) we have defined the new set of 
parameters 



J, = 



(l + a)N c -l 



(J-l) 



b = 



D, 



1 + a 
ea 

1 + a 



(J - J.) 



(1 + ^) 2 
(l + a)(J-l) 



(10) 

(11) 
(12) 

(13) 



We remind in passing that the same equation ([9]) has 
been derived by Willemsen et al. [3(J [H| to describe po- 
larization switches in VCSELs (see also Ref. [35] for a 
similar reduction) . The starting point of their derivation 



is the San Miguel- Feng-Moloney model [36(. The physi- 
cal meaning of the variable (f> is different from here as it 
represents the polarization angle of emitted light. This 
supports the above claim that, upon a suitable reinterpre- 
tation of variables and parameters, many of the results 
presented henceforth may apply also to the dynamics of 
VCSELs. 

In absence of modulation (5 J — 0), Eq. © is bistable 
in an interval of current values where it admits two stable 
stationary solutions (f>± and an unstable one 4>q (double- 
well) . This regime correspond to the bistability region of 
model dl]). Notice that for J = J s , b = the hopping 
between the two modes occurs at the same rate. The 
above definitions allows an estimate of relaxation time Tr 
defined above. This is is the inverse of the curvature of 
the potential in (f>Q. For Jq = J s this is straightforwardly 
evaluated to be 



Tr 



(i + <0 



(14) 



S( J s - 1) 

For the two parameter sets given in Table I one finds 



6 



Tr = 210 Tr = 77.0, respectively. These are the values 
cmploied to draw the leftmost vertical lines in Figs. [TJ31 
The effect of a time-dependent current is to make the 
coefficients a, b and D<p fluctuating. It can be shown [2l[ 
that the effect on can be recasted as a renormalization 
of the intensity of the spontaneous-emission noise. How- 
ever, for the parameters employed in the present work it 
turns out that this correction is pretty small and will be 
neglected henceforth by simply considering as con- 
stant [38]. For simplicity, we also disregard the depen- 
dence of Dc/, on SJ in the drift term of Eq. ([9]). Under 
those further simplifications the Langevin equation can 
be rewritten as 



</> = -U'(<f>)-V'((t>)5J + y/2D^i 



(15) 



where we have express the force term as derivatives of 
the "potentials" 



U{<j>) = 

v{4>) = 



■ COS 1 



S(Jo - 1) 
" 16(1 + a) 
-D<f, In sin2</> 
6 

: r COS 4<» — 

16(1 + 0-) v 



£Q-(Jq - J s ) 

4(1 + a) 



cos 20 



so 



4(1 + a) 



(16) 

cos 20. (17) 



Langevin equations of the form (| 1 5|) with (H} have been 
thoroughly studied in the literature (see e.g. [7l. IqL ITol. 
ITTj ] and references therein) as prototypical examples of 
the phenomenon of activated escape over a fluctuating 
barrier. In view of their non-Markovian nature, their 
full analytical solution for arbitrary t is not generally 
feasible. Several approximate results can be provided in 
some limits. 

For an arbitrary choice of the parameters, V has a 
different symmetry with respect to U meaning that the 
effective amplitude of multiplicative noise is different 
within the two potential wells. If this difference is large 
enough, current fluctuation will remove the degeneracy 
between the two stationary solutions. This is best seen 
by computing the istantaneous potential barriers AU±(t) 
close to the symmetry point Jq = J s . For weak noise 
and S J <C (J s — 1), they are given to first-order in SJ(t) 
by 



AU±(t) 



(1 + *) 



(Js - 1) + 



S±2ea 
8(1 +aj 



6J(t). (18) 



Obviously, this last expression makes sense only when 
the fluctuating term is sub-threshold i.e. whenever the 
system is bistable. In the case of periodic modulation, 
formula (| 1 8|) allows estimating the range of amplitude 
values for a sub-threshold driving 



A < 



SjJs - 1) 
S±2ea 



(19) 



Using this condition, along with the parameter values at 
hand, we deduce that the cases displayed in lower panels 
of Figs. [3] and U correspond to superthreshold driving. 



However, while the minima are much more pronounced 
than in the other panels, there is no qualitative difference 
in the system response. In the case of stochastic modu- 
lation, the same remark applies in a probabilistic sense 
for the last panels of Figs. Q] and [2j 

Altogether, the mode switching can be seen as an acti- 
vated escape over fluctuating barriers given by Eq. (|18[) . 
The statistical properties of the latter process is con- 
trolled by the current fluctuations. We now discuss the 
properties of various regimes. For simplicity, we refer 
to the case of stochastic modulations. Most of the re- 
marks and formulas reported in the following Subsection 
should apply also to the periodic case by replacing r and 
Dj with 2ir/£l and A 2 /2 whenever appropriate. 



A. Fast barrier fluctuations: t < Tr -C T± 

As we already pointed out, in this regime the reduc- 
tion to Eq. (fT5)) is not justified. We may thus only expect 
some qualitative insight on the behaviour of the rate- 
equations. From a mathematical point of view, some 
analytical approximations for equations like (|15p are fea- 
sible in this limit (see e.g. Ref. Q for the stochastic 
case). For our purposes, it is sufficient to note that in 
this regime the effect of 6 J is hardly detected for both 
types of driving (see again Figs. [TJ4]) . Note also that 
working at Dj fixed means that for r — > the fluctua- 
tion become negligible. 



B. Resonant activation: Tr < t <C T± 

If Tr < t we are in the colored noise case. The prob- 
lem is amenable of a kinetic description which amounts 
to neglect intrawell motion and reduce to a rate model 
describing the statistical transitions in terms of transition 
rates. If we consider r as a time scale of the external driv- 
ing we can follow the terminology of Ref. [39( and refer 
to this situation as the "semiadiabatic" limit of Eq. (|15p . 

In this regime, the residence time is basically the short- 
est escape time, which in turn correspond to the lowest 
value of the barrier (the noise is approximatively con- 
stant in the current range considered henceforth). For 
the case of interest, 5 < 2ea we can use (fT5|) to infer 
that the minimal values of AU± should be attained for 
SJ cx +\ADj respectively. This yields 



T± ~ T s exp 



K- 



2ea±5 



1 



Dj 



(20) 



where K is a suitable numerical constant. Notice that 5 
controls the asymmetry level: if 8 <C 2ea the two resi- 
dence times decrease at approximatively the same rate. 
This prediction is verified in the simulations and also in 
the experiment [2lj |. 

As a further argument in support of the above rea- 
soning, we also evaluated the probability distributions 
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of the residence times obtained from the simulation of 
the rate equations. In Fig. [5J we show two representa- 
tive cumulative distributions. The data are well fitted 
by a Poissonian P(T) = 1 - exp(-T/T±) for both the 
stochastic and periodic modulation cases. This confirms 
that hopping occours preferentially when a given (mini- 
mal) barrier occurs. 





FIG. 5: (Color online) Cumulative distributions of the resi- 
dence times in the resonant activation region, parameter set 
with c = 1.3 (see text and Table I). Left panel: stochastic 
modulation with Dj = 5xl(T 4 , r = 1.638 x 10 3 . Right panel: 
periodic modulation with A = 0.03 and period 1.286 x 10 4 . 
We report only the histograms for the times whose averages 
are denoted by T+ in the text. Solid line is the cumulative 
Poissonian distribution with the same average. 



C. Slow barrier, frequent hops: Tr <C T± <C r 

This corresponds to the adiabatic limit in which the 
time scale of the external driving is slower than the intrin- 
sic dynamics of the system [39( . To a first approximation 
we can here treat current fluctuations in a parametric 
way. Correction terms may be evaluated by means of 
a suitable perturbation expansion in the small parame- 
ter 1/t If S J is small enough for the expression 
(ITS)) to make sense, the escape time can be estimated as 
the average of escape times over the distribution of bar- 
rier fluctuations, i.e. (T±)sj- For the case of Eq. ([J|, 
the variable 5 J is Gaussian and we can use the identity 
(expf3z) = exp(/3 2 (z 2 )/2) to obtain [TO] 



T± ~ T s exp 



2(S±2ea) 2 



Dj 



(21) 



This reasoning implies that for large r the residence times 
should approach two different constant values. A closer 
inspection of the graphs (in linear scale) reveals that this 
is not fully compatible with the data of FigfT] even for 
the smallest value of Dj. In several cases, T± continue to 
increase with r and no convincing evidence of saturation 
is observed. We note that the same type of behaviour 
was already observed in the analog simulations data of 



Ref. [ll|. There, an increase of hopping times duration at 
large r was found. The Authors of Ref. [ll[ explained this 
as an effect of a too large value of the noise fluctuation 
forcing the system to jump roughly every r. We argue 
that the same explanation holds for our case. This is also 
consistent with the fact that the exponential factors in 
Eq. (|21|) evaluated with the simulation parameters turn 
out to be much larger than unity. 



V. CONCLUSIONS 

In this paper, we have explored numerically and ana- 
lytically the effects of external current fluctuations on the 
mode- hopping dynamics in a model of a bistable semicon- 
ductor laser. To the best of our knowledge, this setup 
provides the first theoretical evidence of resonant activa- 
tion in a laser system. As the phenomenon has hardly 
received any experimental confirmation in optics, we be- 
lieve that our study may open the way to future research 
in this subfield. 

The model we investigated is based on a rate-equation 
description, where the bias current enters parametrically 
into the evolution of the modal amplitudes. We consid- 
ered, two kinds of current flutuations, namely, a stochas- 
tic process ruled by an Orstein-Uhlenbeck statistics, and 
a coherent, sinusoidal modulation. These choices are mo- 
tivated by the aim of proposing a suitable setup for an 
experimental verification of our results. Upon varying the 
characteristic time-scale of the imposed fluctuations, we 
have shown that the residence times attain a minimum 
for a well-defined value, which is the typical signature 
of resonant activation. The magnitude of the effect can 
be different depending on the parameters of the model. 
Moreover, the response of the system appears very much 
similar for both periodic and random modulations. 

The reduction of the rate equations to a one- 
dimensional Langevin equation allowed us to recast the 
problem as an activated escape over a fluctuating bar- 
rier. To first approximation, the fluctuating barrier (mul- 
tiplicative term) is mainly controlled by current modu- 
lations while the spontaneous noise act as an additive 
source. This simplified description has allowed us to 
draw some predictions (e.g. the dependence of residence 
times on noise strength) and to better understand the 
role of the physical parameters. Given the generality of 
the description, our results should apply to a broad class 
of multimode lasers, including both Edge Emitting and 
Vertical Cavity Lasers. 

From an experimental point of view, driving the laser 
in a orders-of-magnitude wide range of time-scales is 
more feasible in the case of a sinusoidal modulation than 
for a colored, high frequency noise. However, given the 
evidence of a resonant activation phenomenon for such 
modulation, our results indicate that it occurs almost for 
the same parameters in the case of colored noise, provided 
that the RMS of the modulations equals the amplitude 
of the added noise. Thus, the phenomenon could be fully 
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exploited along those lines. Since the reported experi- characterization in optical systems that allows for both 
mental evidences of the phenomenon are so far scarce, very precise measurements and careful control of param- 
we hope that the present work could suggest a detailed eters. 
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